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Abstract: Recently, LHCb experiment announced a 3.5a evidence for direct CP violation in D° decay by 
looking at the difference between Acp(D° — > K + K~) and Acp(D° — ¥ tt + -k~). This is the first evidence of 
CP violation in charm system, which may indicate new physics beyond the Standard Model. Motivated by 
this measurement, we review rare processes in charmonium decay, especially, the weak decay, C or P violated 
decay and lepton flavor violated decays. In case the new physics appears in charm sector, these rare decays 
of charmonium states will provide opportunity to search for significant contributions from physics beyond the 
Standard Model. With huge J/ip and tp(2S) samples in BESIII experiment, the rare decays may be feasible. 
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1 Introduction 



At present, two general trends can be distin- 
guished in accelerator particle physics. On one hand, 
very high energy accelerators, for example, the LHC, 
provide the ability to explore physics at high energy 
frontier. On the other hand, smaller experiments at 
lower energies but with very high intensities and low 
backgrounds, for example, the B factories and BE- 
SIII , provide the capabilities for performing precise 
tests and accurate determinations of many parame- 
ters of the Standard Model (SM). Moreover, the close 
scrutiny of rare processes may illuminate new physics 
in a complementary fashion to high-energy colliders. 

On November 14, 2011, LHCb experiment an- 
nounced 3.5a evidence for direct CP violation in 
D° decay. They looked at the difference between 
A C p(D° K+K-) and A CP (D° -> tt+tt-), and 
found that the AA CP = [-0.82 ± 0.21 (stat.) ± 
0.11(sys.)]% Q. Since the SM predicted that the 
CP asymmetry in charm sector should be less that 
0.1% [2j, this is the first indication of new physics 
at LHC. If the new physics could appear in the up 
quarks (charm and top), it urgently needs confirma- 
tion. One possible way is to probe the charmonium 
rare decays with huge J/ip and ip(2S) data sample. 
The BESIII experiment can provide important tests 
of the SM, with the accompanying possibility for un- 



covering new physics induced deviations in charm sec- 
tor, especially, the weak charmonium decay, invisible 
decay, lepton flavor violated decays and so on. In this 
paper, we review the possibility of such kinds of rare 
charmonium processes. 



2 Weak decays of charmonium 



The low-lying charmonium states, i.e. those be- 
low the open-charm threshold, usually decay through 
intermediate photons or gluons produced by the an- 
nihilation of the parent cc quark pair. These OZI- 
violating but flavor-conserving decays result in nar- 
row natural widths of the J/ip and ip(2S) states. In 
the SM framework, flavor-changing weak decays of 
these states are also possible, although they are ex- 
pected to have rather low branching fractions. The 
huge J I ifi data samples at BESIII will provide oppor- 
tunities to search for such rare decay processes, which 
in some cases may be detectable, even at SM levels. 
The observation of an anomalous production rate for 
single charmed mesons in J/ip or ip(2S) decays at BE- 
SIII would be a hint of possible new physics, either in 
underlying continuum processes via flavor-changing- 
ncutral-currents 0-0] or in the decays of the ip res- 
onances due to unexpected effects of quark dynam- 
ics [13, Eli- 
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2.1 Semileptonic decays of charmonium 

The inclusive branching fraction for J / tp weak de- 
cays via a single quark (either the c or the c) had 
been estimated to be (2 ~ 4) x 10~ 8 by simply using 
the D° lifetime Such a small branching fraction 
makes the observation of weak decays of the J/ip or 
ip(2S) quite challenging, despite the expected clean- 
liness of the events. However, BEPC-II, running at 
designed luminosity, will produce of order 10 10 J/ip 
events per year of data taking, leading to = 400 weak 
decays for the predicted SM branching fraction. The 
semi-leptonic decay of a cc (1~~) vector charmonium 
state below the open-charm threshold is induced by 
the weak quark- level transition c — > qW* , where W* 
is a virtual intermediate boson. Hence, the accessible 
exclusive semi-leptonic channels are: 



ip(nS) — > D q lv, 



ip(nS) - 



■ d;iu, 



(i) 



(2) 



where n = 1 or 2, and q can be either a d- 
or s-quark, which corresponds to a Z3 ± (Cabibbo- 
suppressed mode) or D s (Cabibbo-allowed mode) me- 
son. Semi-leptonic weak decays of the J/ip will of- 
fer several advantages over the purely hadronic ones 
from both the experimental and theoretical points 
of view: the prompt charged lepton I = e,/j, can 
be used to tag the events, removing a large frac- 
tion of conventional ip(nS) hadronic decays. In addi- 



tion, the missing energy due to the escaping neutrino 
can be also exploited to remove backgrounds. The 
identification of the charm meson in the final state 
would provide an unambiguous signature of the semi- 
leptonic weak decays of ijj(nS). Meanwhile, decays of 
the excited mesons D* and D r± produced in reac- 
tion would provide useful additional experimental 
handles. In the lab system, the detectable photons 
from the D** — > Dfj, radiative transition are in the 
90~200 MeV energy interval. These, and the soft 
pion produced from D* ± in D* ± — )• D°tt ± decay, can 
provide powerful constraints to help identify a D 3 or 
D° meson produced in the weak decay of a charmo- 
nium state. 

Many theoretical calculations have been done 
based on various QCD frameworks, including, e.g., 
heavy quark spin symmetry (HQSS) [4j, QCD sum 
rules Q and covariant light- front quark model Q. 
The predicted branching fractions are in the range of 
10- 9 -10- 10 and 10- 10 -10- n for J/tp -> D [ * ] lv and 
J /ip —¥ D^*'lf, respectively. The branching fractions 
for the sum of dominant semi-leptonic decay modes 
can arrive at 5 x 10~ 9 Q, which are hopefully to be 
detected at BESIII with 10xl0 9 J/i/j events accumu- 
lated each year. The overall theoretical uncertainty 
is up to ~ 50% due to the treatment of the nonper- 
turbative QCD dynamics. In Table [TJ the theoretical 
predictions for J/?p semi-leptonic decays are summa- 
rized. 



Table 1. The predicted branching fractions (in units of 10 -10 ) for J/ip semi-leptonic weak decays. The 
transition form factors for J/ip—¥ D^* d are estimated based on the ISGW model in Ref. Q], the factorization 
scheme in Ref. [f| , the QCD sum rules in Ref. 0] and the covariant light-front quark model in Ref. Q . Here, 



BK(J/il> ->■ D*-l+u) = B1Z{J/iP D* s - d e+u e )+BTZ(J/^ ■ 
the theoretical predictions. 



Only the central values are cited from 



decay mode 


Ref. [4] 


Ref. [5] 


Ref. [6] 


Ref. [7] 


J/ifj^Djl+v + c.c. 


26 


20 


3.5 


11.2 


J/V>-» D* a -l + v + c.c. 


42 


32 


11 


34.7 


J/V>-» D~l + v + c.c. 


1.4 


1.2 


0.14 


1.1 


J/ijj^D*-l + v + c.c. 


2.3 


2.0 


0.73 


3.4 



It is interesting to note that the ratio of Cabibbo- 
allowed decay to Cabibbo-suppressed decay can be 
obtained cleanly since many theoretical uncertain- 
ties cancel out in the ratio. The ratios R s /d — 

bk(J/^d:-i + v) shou i d be 
BTz(j/i>^D*-i+v) bmjulu un- 
equal to j^jT — 18.4 under the SU(3) flavor sym- 
metry limit [J], where V ca = 0.973 and V cd = 0.2271 
denote the relevant Cabibbo-Kobayashi-Mashkawa 



BTZ(.J/^^D-l + u) , R » 

BTl(.J/i,^D-l + v) cU1U ^a/d 



(CKM) mixing matrix elements. In Ref. [6j, the cal- 
culations based on the QCD sum rules show R s/d = 



24.7 and R 



Id 



15.1, which implies a large effect 



of the SU(3) symmetry breaking. Under the as- 
sumption of heavy-quark spin s ymm etry (HQSS) and 
the non-recoil approximation [12l 

,-, BK(J/i)^D*-l+is) 



13| , the ratio of 



W/P 



is predicted to be 1.6 



BKiJ/ip-tD'l+v) 

while it is calculated to be 3.1 based on the QCD 



No. X 



Li Hai-Bo, Zhu Shi-Hai: Mini-review of rare charmonium decays at BESIII 



3 



sum rules in Ref. which indicates that high order 
corrections are important. 



Following the first proposal by Zhang 14j , BESII 
has measured the J/ip weak decays to be BlZ(J/ip^>- 
D~e+v + c.c.) < 3.6 x 1CT 5 and BTZ{J/ip -> D~e+u + 
c.c.) < 1.2 x 10~ 5 at 90% confidence level based on 
5.8 x 10 8 J/ip decay events The current exper- 

imental limits are much higher than the predicted 
values. It will be important to search for J/ip weak 
decays at the BESIII experiment, which will be dis- 
cussed in detail in Section 2.3. 

The flavor-changing-neutral-current (FCNC) in- 
duced semi-leptonic decays of J/ip —> L>^* )0 Z + Z~ are 
studied in Ref. § based on the QCD sum rules, and 
the decay rates for J/ip -> D°l+l' and J/ip -> D*°l + l' 
are predicted to be 10~ 14 and 10~ 13 , which are much 
below the BESII upper limit: BTZ(J/ip -> D°e+e-) < 

1.1 x 10~ 5 [15]. But new physics in the loop may 
enhance the decay rate considerably as discussed in 
Ref 14 1, further experimental search at BESIII will 
be interesting. 

2.2 Two-body weak hadronic decays of char- 
monium 

Non-leptonic two-body weak decays of the J/ip 
and ip(2S) are studied based on the covariant light- 
front quark model 0] , the factorization approach Q , 
and HQSS 0, Q in the context of the factorization 
scheme for both the Cabibbo-allowed (c — > s) and 
Cabibbo-suppressed (c — > d) quark-level transitions. 



In the literature, the expressions for branching frac- 
tions for J/ip -> PP/PV/VV decays (where P and V 
represent pseudoscalar and vector mesons) had been 
given. For J/ip — > PP mode, the dominant decay is 
J /ip — > D~n + + c.c. (including charge conjugate pro- 
cess) , which is Cabibbo and color allowed. Table [2] 
shows the decay rate for the J/ip — > PP modes from 
various theoretical predictions EBB El. 

Assuming factorization as suggested by 
Bjorken [17], it is interesting to construct the ra- 
tio of Cabibbo suppressed to Cabibbo allowed decay 
modes: 



BK{J/iP^D s K) \V ua f R 



BK(J/^j^D 3 n) \V ud f„ 



' 0.081, 



(3) 



where Jk and f n are the decay constants for K and 
7r mesons. The ratio was predicted to be 0.08 based 
on the factorization approach Q. 

For J/ip — > PV mode, the Cabibbo favored and 
color allowed decay J/ip —> Dfp~ +c.c. is dominant 
and predicted to be (12 ~ 50) x 10~ 10 as shown in Ta- 
ble^ which can be accessible in BESIII experiment. 
Moreover, the ratio, ^ J/ ^^ sP \ was predicted to 
be 6.3 in [9j and 4.2 in 16j, which can be tested by 



BESIII experiment. 

For J/ip -> PV mode as listed in Table IH the 
dominant decay mode is J/ip —¥ D* a + p~ +c.c, which 
was predicted to be 5.3x 10~ 9 based on the the fac- 
torization approach Q . This decay mode is the most 
promising mode to be measured at BESIII. 



Table 2. The theoretical predicted branching fractions (in unit of 10 -10 ) for ip^> PP . The transition mode, 
AC = AS = +1, corresponds to the Cabibbo-allowed decay modes, while AC = +1, AS = corresponds to 
the Cabibbo-suppressed decay modes. Only the central values are cited from the theoretical predictions. 



Transition mode 


Decay mode 


Ref. [5] 


Ref. ]7] 


Ref. ]9] 


Ref. [16] 


AC = A5 = +l 














ip — > D~ir + +c.c. 


7.36 


2.5 


2.0 


8.74 




ip^D°K° + c.c. 


1.39 


0.5 


0.36 


2.80 


AC = +1, AS = 














ip — > DfK~ +c.c. 


0.53 




0.16 


0.55 




ip — > D +/ n~ -\-c.c. 


0.29 




0.08 


0.55 




ip -> D°rj + c.c. 


0.070 






0.016 




ip->D°ri' +c.c. 


0.004 






0.003 




1p^D°TT°+C.C. 


0.024 






0.055 



The J/ip semi-leptonic decay modes can be re- 
lated to the two-body hadronic decay modes by ap- 
plying both the spin symmetry and the non-recoil ap- 
proximation to the semi-leptonic decay rates Q. For 



J/ip D+(D* s + )ir~ decay modes, q 2 = (py, -p D ) 2 = 
m 2 (here and p D are the four momenta of the ini- 
tial and final state heavy mesons) and, assuming fac- 
torization as suggested by Bjorken [17| for B decays, 
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and in the non-recoil approximation for the hadronic 
transition amplitudes [18], one can give the relation 
between relative branching ratios: 

BR{J/i>^>D*+%-) 



BR(J/ip- 

dr(j/ip- 



■Dl+l-is)/dq 2 



dT{J/iP^D+l-v)/dq 2 



(4) 



which was predicted to be 7.5 based on the the fac- 
torization approach 9]. If a p is substituted for the 
7r one gets r = 4.2 in 9j. In this way, the estimated 
branching ratios in Table [5] for if)(nS) — > PP chan- 
nels can be related to if)(nS) — > VP channels with the 
pseudoscalar charm mesons replaced by vector charm 
mesons. 



Table 3. The theoretical predicted branching fractions (in unit of 10 -10 ) for J/if)—¥ PV . The transition mode, 
AC — AS — + 1, corresponds to the Cabibbo-allowed decay modes, while AC = +1, AS = corresponds to 
the Cabibbo-suppressed decay modes. Only the central values are cited from the theoretical predictions. 



Transition mode 


Decay mode 


Ref. [5] 


Ref. [7] 


Ref. [9] 


Ref. [16] 


AC = AS=+1 














ip — > Df p~ + c.c. 


50.5 


28.0 


12.6 


36.30 




iP^D K*° + c.c. 


8.12 


5.5 


1.54 


10.27 


AC = +1, A5 = 














ip^D+K*~ +c.c. 


2.79 




0.82 


2.12 




if) — ► D + p~ +c.c. 


2.13 




0.42 


2.20 




ip->D°p° + c.c. 


0.18 






0.22 




tp -> D°LU + C.C. 


0.16 






0.18 




D°cf> + c.c. 


0.41 






0.65 



Table 4. The theoretical predicted branching 
fractions (in unit of 10" 10 ) for J /if) ->■ VV 
. Only the central values are cited from the 
theoretical predictions. 



Decay mode 


Ref. [9] 


J/ip^-D*+p-+c.c. 


52.6 


J/ip^D*+K*- +c.c. 


2.6 


J/if^D'+p- +c.c. 


2.8 


J/ip->D*+K*-+c.c. 


9.6 



All of the above estimates show an overall en- 
hancement of the final-state vector charm mesons 
with respect to the pseudoscalar ones. This suggests 
the use of D* or D* ± as signals in searches for weak 
decays of the J /if in non-leptonic decay channels as 
well. According to the predicted results in Ref. 
the branching fraction for inclusive weak hadronic de- 
cays of J /if) can be as large as 1.3x 10~ 8 , which is in 
remarkable agreement with the naive estimation in 
Ref.fl. 

In the SM, FCNC induced J /if weak hadronic de- 
cays are predicted to be unobservably small [Tol ] and, 
thus, any observation of such decay would provide a 
signal for new physics. In Ref. [lfj the predictions 
of various models, such as TopColor models, mini- 
mal supersymmetric standard model (MSSM) with 



R-parity violation and a general two-Higgs-doublct 
model, are discussed. These authors found that the 
branching fraction for J /if) — > D/DX U , which is me- 
diated by the weak c — > u transition, could be as large 
as 10 -6 ~ 10 -5 in some new physics scenarios. 

At BESIII it will be difficult to isolate pure, c — > u 
mediated, hadronic J /if) — > Z? w /D M X U decays. On 
the other hand the decays J/tf) D 0M / L> ow 1+ 1~ 
(l — e, p) and J/ip-> D 01 ^ /D° M j decays, which are 
also dominated by FCNC processes, would be quite 
distinct. 

In addition to the J /if) and if>{2S) decays into 
charm meson, the weak decays of ip(2S) into final 
states involved A c baryon can also be searched for at 
BESIII. These decay modes include if(2S) -> A+S" + 
c.c, A+S°7r" +c.c, A + S _ e~P+c.c and A+A7r~+c.c. 
In ref. [19], the decay of ^(25) ->• A+S~ +c.c. are es- 
timated in the framework of the SM, and they found 
that the decay rate could reach 10~ 10 . 

2.3 Searches at BESIII 

At BESIII, assuming a 10 10 J /if) event sample, 
the total semi-leptonic decay rate could be 5 x 10~ 9 , 
so about 50 semi-leptonic decay events of the type 
J/tf) — > D a (D* a )lv can be accumulated by BESII per 
year. The following event selection criteria would be 
useful for searching such exclusive semi-leptonic chan- 



No. X 



Li Hai-Bo, Zhu Shi-Hai: Mini-review of rare charmonium decays at BESIII 



5 



nels : 

• The prompt charged lepton can be used to tag 
the weak decay: in order to suppress cascade 
decay backgrounds from J /tp strong decays, the 
tagging lepton momentum could be required to 
be between 0.5 GeV and 1.0 GeV, close to the 
upper kinematic limit for the decay under con- 
sideration. High quality lepton discrimination 
from charged pions or kaons is needed for the 
measurement. 

• The missing mass of the reconstructed candi- 
dates must be consistent with the (nearly) zero 
mass of the undetected neutrino. 

• The reconstruction of a D a or D meson would 
provide an unambiguous signature for a weak 
decay of a below-open-charm threshold ijj(nS). 
Good invariant mass resolution of the D s decay 
products will be important for reducing combi- 
natorial backgrounds. 

• Soft photons in the energy interval 
(90-200) McV from the D* g ->■ jD 3 transition 
and soft charged pions from D f± — > 7r ± D decay 
can provide further suppression of combinato- 
rial backgrounds: the additional constraint of 
an intermediate D* state would reconfirm the 
D s signal. 

In general, exclusive hadronic decays are probably 
too tiny to look for in any specific fully reconstructed 
decay channel. Therefore, it seems that an inclusive 
search for J/ip — y D* + X at BESIII may be more 
fruitful. The 7 from the decay of a D* meson should 
be useful as a kinematic constraint to clean up any 
D s meson signal, as discussed in 

3 Search for the invisible decays of 
quarkonium 

Invisible decays of quarkonium states offer a win- 
dow to look for possible new physics beyond the 
SM 2(3, 21 1 . The reason is that other than neutrinos, 
the SM includes no other invisible final particles that 
these states can decay into. BESII explored such a 
window by establishing the first experimental limits 
on invisible decays of the 77 and rj', which comple- 
mented the limit of 2.7 x 10 -7 that was previously 
established for invisible decays of the ir° [221 ] . 

Some theories of beyond the SM physics predict 
new particles with masses that are accessible at BE- 
SIII, such as the light dark matter (LDM) particles 



discussed in Ref. 23[. These can have the right relic 
abundance to constitute the nonbaryonic dark matter 
of the universe, if they are coupled to the SM via a 
new light gauge boson U 2J| , or exchanges of heavy 
fermions. A light neutralino with a coupling to the 
SM that is mediated by a light scalar singlet in the 
next-to- minimal supersymmetric standard model has 
also been considered [251 ] . 

These considerations have received a boost in in- 
terest by the recent observation of a bright 511 keV 
7-ray line from the galactic bulge reported by the 
SPI spectrometer on the INTEGRAL satellite Q. 
The corresponding galactic positron flux, as well as 
the smooth symmetric morphology of the 511 keV 
7 emission, could be interpreted as originating from 
the annihilation of LDM particles into e + e~ pairs 



[23L 1271 ] . It would be very interesting to see evi- 
dence for such light invisible particles in collider ex- 
periments. CLEO gave an upper bound on T(1S) — > 
7 + invisible, which is sensitive to dark matter can- 
didates lighter than about 3 GeV/c 2 [2^], and also 
provided an upper limit on the axial coupling of any 
new U boson to the b quark. It is important, in ad- 
dition, to search for the invisible decays of othr light 
quarkonium states (qq, q — u,d, or s quark), since 
these can be used to constrain the masses of LDM 
particles and the couplings of a U boson to the light 
quarks [21 1. 

It has been shown that measurements of the J / ip 
invisible decay widths can be used to constrain new 
physics models It is straightforward for one to 

calculate the SM ratio of branching fractions for J / ip 
invisible decays and its measured branching fraction 
for decays into electron-positron pairs 29(. Within 
the SM, the invisible mode consists solely of decays 
into the three types of neutrino-antineutrino pairs. 
Neglecting polarization effects and taking into ac- 
count e + e~ production through a photon only, one 



gets [29|: 

r(j/^^e+e-; 



27G 2 Mj /ii 
2567r 2 a 2 
= 4.54 xlO" 7 



(5) 



where G and a are the Fermi and fine-structure con- 
stants, respectively, and Mj/^ is the J/ip mass. The 
uncertainty of the above relation is about 2—3% and 
comes mainly from the corrections to the J/ip wave 
function, e + e~ production via the Z boson and elec- 
troweak radiative corrections [29| . 

In BES experiment, one can tag the charmonium 
states that decay invisibly by looking for a particu- 
lar cascade transition, such as ip(2S) —> ir + ir~J/tjj, 
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ip{2S) — s> 7% c and so on, where the soft 7r + 7r~ pair or 
the monoenergetic radiative 7 serves as a tag for the 
invisibly decaying J/ip or \c state. The BESII exper- 
iment performed the first search for invisible decays 
of the J/ip using tp' — > Tr + ir~J/ip events detected in a 
sample of 14.0 million tp' decays. The upper limit on 
the ratio B( ^f ^ inv j sib | e) at the 90% confidence level 
is 1.2 x 10 -2 [30]. This measurement improves by a 
factor of 3.5 the bound on the product of the coupling 
of the U boson to the c quark and LDM particles as 
described in Eqs. (25) and (26) of Ref. 31] . One now 



has, for a Majorana LDM particle x as m Eq. (26) 
of Ref. [il|, a limit of \c x f cV \ < 8.5 x 10~ 3 , which is 
almost a factor of 2 stronger than the corresponding 
limit |c x /;,v| < 1.4 x 10~ 2 derived from the invisible 
decays of the T(1S) as described in Eq. (106) in 
Ref. [32j, where c x and f cV (fbv) denote the U boson 
couplings to the LDM particle x an d c (6) quark. We 
expect a more precise measurement can be obtained 
in the future BES-III experiment. A list of poten- 
tially useful decay chains is provided in Table [5] in 
the BESIII experiment. 



Table 5. ip(/2S) and J/ip decay modes that can be used to search for invisible decays of the J/ip, XcO, Xm, 
Xc2, ti c (1S) and t) c (2S). The branching fractions are taken from the PDG [3^]. For each mode, a "tagging 
topology" is given, which is the set of visible particles that are seen within the detector's acceptance. In 
each case the tagging topology has well defined kinematics. The number of events are the expected event 
yield in 3 billion ip(2S) (10 billion J/ip) data set, in which we did not consider the decay probabilities of the 
tagging particles. 



■p{2S) 


Branching 


Number of events 


Invisible 


Tagging 


decay mode 


fraction (10~ 2 ) 


expected 


decay mode 


topology 


V>(2S) ^tt+tt" J/ip 


31.7±1.1 


9.3 x10 s 


J/ip—t invisible 


7T~^7T — 


V>(2S)^7r 7r°J/»/> 


18.6 ±0.8 


5.6 x10 s 


J/ip—t invisible 




ip(2S)-*T)J/ip 


3.08±0.17 


9.3 xlO 7 


J/ip—t invisible 


V 


ip(2S)^n°J/iP 


0.123±0.018 


3.7 x 10 6 


J/ip—t invisible 


7T 


V>(2S)->-7XcO 


9.0±0.4 


2.7 x 10 s 


XcO - ► invisible 


7 


ip(2S)-+-rXcl 


8.7±0.5 


2.6 x10 s 


Xci invisible 


7 


V>(2S)-> 7 Xc2 


8.2±0.3 


2.5 x10 s 


Xc2 invisible 


7 


■p(2S)^rnc(is) 


0.26 ±0.04 


7.8 x 10 6 


r) c (lS) — > invisible 


7 


J/^^7Vc(lS) 


1.3 ±0.4 


1.3 x10 s 


r; c (lS) — > invisible 


7 



It is also interesting to search for invisible de- 
cays of the 77, 77', p, lo and <p light mesons. Within 
the SM, the decays of 77(77') — > vv are tiny due to 
the helicity suppression. If the Z° couples to a 
massive neutrino with the standard weak interaction 
strength, the branching ratios (BR) for n° — > v T v T 
and 77 — > v T i> T have maximum values of 5.0 x 10 -10 and 
1.3xl0 -11 [3J,|35(, respectively, at the v T mass upper 
limit of m VT = 18.2 MeV/c 2 33j. Any enhanced signal 
of invisible decay may indicate New Physics. Possible 
77/7/ — ^invisible decay products could be LDM parti- 
cles or light neutralinos. These LDM particles may 
have an adequate relic density to account for the non- 
baryonic mass of the universe. 

In the BES experiment, the two-body decay 
modes J/ip — > ^77 or (prj' can be selected using only 
the very clean and distinct (j> — > K + K~ decays, which 
then tag the presence of an 77 or 77' meson that has 
decayed into an invisible final state. For J/ip <prj' , 
the missing momentum, P misa — |P miss |, is a power- 
ful discriminating variable to separate signal events 
from possible backgrounds in which the missing side 



is not from an 77 (77'). Here, P misa = —P^,. In addition, 
the regions of the detector where the 77 and 77' decay 
products are expected to go are easily defined thanks 
to the strong boost from the J/ip decay. 

The branching fraction of 77(77') — > 77 was also 
determined in J/ip — > 4>rj{rj') decays, in order to pro- 
vide the ratio 6(77(77') — > invisible) to £> (77(77') — > 77). 



The advantage of measuring 



6(77(77') — > invisible) . 



8(77(77') -> 77) 
that the uncertainties due to the total number of 

J/ip events, tracking efficiency, PID, the number of 
the charged tracks, the cut on M(KK), and residual 
noise in the BSC all cancel. 

Based on 58 x 10 6 J/ip events in the BESII ex- 
periment, the upper limit on the ratio of the 6(7 7 — > 
invisible) to B(r\ — > 77) was determined to be |36| : 



6(77(77') — > invisible) 
6(77(77') -> 77) 



<1.65xl0" 3 (6.69xl0" 



(6) 



Table [6] lists the possible two-body J/ip decay 
modes that can be used to study the invisible decays 
of the 7/, 77', p, u) and <p mesons at BESIII. 
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Table 6. J/ip decay modes that can be used to study invisible decays of rj, rj ', p, ui and (p mesons. The 
branching fractions are from the PDG 33,!]. For each mode, a "tagging topology" is given, which is the set 
of visible tracks in the detector's acceptance. In each case the tagging topology has well defined kinematics. 
The produced number of events are the expected events in 10 billion J/ip event data set at BESIII, with the 
decay probabilites of the tagging particles included. 





Branching 


Invisible 


Tagging 


Number of events 


decay mode 


fraction (1CT 4 ) 


decay mode 


topology 


/10 billion J/4> decays 


J /if) — y (jyq 


6.5±0.7 


Tj — y invisible 


<f>^K+K- 


(31.4±3.4) X 10 5 




6.5±0.7 


<f> —y invisible 


Tj— yyy 


(25.7±2.8) x 10 5 


J/ip — > (fir)' 


3.3±0.4 


rf — y invisible 


<j)^K+K- 


(16.2±1.9) x 10 5 




3.3±0.4 


4> — ¥ invisible 


rf -y 7p° 


(9.6±1.2) x 10 5 


J/lf) — > LIT] 


15.8±1.6 


Tj — > invisible 




(13.9±1.4) x 10 6 




15.8±1.6 


ui —y invisible 


rj — y 77 


(6.2 ±0.6) x 10 6 


J/ip -y Lurf 


1.67±0.25 


rj' —y invisible 


UJ — y 7T° 


(1.5±0.2) x 10 6 




1.67±0.25 


to — y invisible 


rf -y 7P° 


(0.7±0.1) x 10 6 


J/ip-y p°V 


1.93±0.23 


•r) — > invisible 


p° —y n+n~ 


(1.9±0.2) x 10 6 




1.93±0.23 


p° — y invisible 


rj— )-77 


(0.8 ±0.09) x 10 e 


J/lp^p°TT° 


56±7 


p° — > invisible 


tt° — y 77 


(55.3±5.8) x 10 6 



Search for C or P violating processes 
in J/ip decays 



With its huge Jji\> and ip(2S) data samples, the 
BESIII experiment will be approaching the statistics 
regime where studies of rare ip decays can provide 
important tests of the SM and possibly uncover devi- 
ations. Among the interesting examples are C, P or 
CP violating processes in J/ifj decays. An example 
of such modes would be ip(nS) — > V°V°, where V^ is 
used to denote J PC = 1~" vector mesons (cf>, u>, p° and 
7) . A distinct signal for this class of event would be 
il>(nS) -> <j}(/} detected in ^(nS) K + K~ K + K~ final 
states. Because of the C violation, ip(nS) V°V° de- 
cays can only occur in the SM via cc annihilation via 
a Z° or W-exchange decays as discussed in Ref. [37 1. 
The rate for this type of weak decay can provide a 
measurement of the charmonium wave function at the 
origin [371 ]. 

In order to make a rough estimate for the rate, 
we first consider just the rate due to the W- 
exchange contribution, which is straightforward to 
compute [38| 



can expect that the SM exclusive BR(tp(nS) - 
rate should be below the level of 10~ 8 and probably 
out of reach of the BESIII experiment. 

Experimentally, there are some possible back- 
grounds that will dilute the signal for J/ip —> <p4> de- 
cays. One major background is J/ip — y "f4><fi, which 
is mainly from J/ip — > jrj^lS), r/ c (lS) — > (p<p. This 
background can be removed by doing a constrained 
kinematic fit. A detailed calculation had been done to 
estimate the background from J/ip —y j(p<p 3jJ • An- 



T(J/iP^ssY 
T(J/ip^ye+t 



m w 



(7) 



where rrij/^, and m w are the masses of J/ip and W bo- 
son, respectively. This leads to BR(J/ip — > ss) weak = 
10 -7 for this weak contribution. To form the (pep final 
state, another ss pair must be produced from the vac- 
uum and these s-quarks have to bind with the out- 
going ss from the cc decay to produce the (p(p final 
state 37]. When this is considered, it seems that one 



other background appears if one studies only 2(KK) 
invariant pair mass distributions. It arises from the C 
and P-conserving reaction J/ip—> <p(KK) S - wav<11 due 
to the fact that the (p mass is only two iS-wave-widths 
away from the KK S'-wave resonance mass, for ex- 
ample, / (980) -> KK. Although it may be difficult 
to subtract in a small statistical sample, one can, in 
principle, remove this kind of background by either 
a spin-parity analysis of the KK pairs in a narrow 
window about (p mass, or by a subtraction normal- 
ized to an observed S'-wave mass peak. To avoid the 
S'-wave contribution, one can reconstruct one (p from 
K + K~ and another <p from the K$K L mode, which 
is not allowed to form an S-wave. It will be easy to 
look for the missing mass of one (p reconstructed from 
K + K~, to see if there is any peak under the (p mass 
region by also requiring K s and K L information in 
the final states. 

It is noted that there is possible continuum back- 
ground produced via a two-photon annihilation pro- 
cess. It is a peaking background that cannot be re- 
moved without considering detailed angular distribu- 
tions in a high statistics sample. It is very hard to 
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deal with this kind of peaking source with a small 
sample of signal events. One way is to use off-peak 
data which are taken below the J / ip mass peak. The 
e + e~ — >77 process has been investigated before [39j |. 
and it has a unique production angle (#*) distribu- 
tion, which is defined as the angle between <p and e~ 
beam direction in the Center-of-Mass (CM) frame. 
The production angle distribution for the two real 
photon annihilation process has the form of 

a{cos9 ) e+e -^ = — — — , (8) 

while, in the process of two virtual photon into V°V° 
pairs, the distribution is (to first order) (4oj |: 

1+COS0* 2 



place indirect limits on CLFV J/ip branching frac- 



cr(cOS#*) e + e -_ >7 * 7 *_ ! .v'0 V o 

where factor k is: 

k = 



k 2 — i 



2m 2 v0 -S 



y/S 2 -4Sr 



(9) 



(10) 



where S is the square of CM energy. In principle, by 
using an angular analysis, one can remove the peak- 
ing background with high statistic data sample. To 
avoid the peaking background from the continuum, 
ip{2S) — > irirJ/ip could be used to study this kind of 
rare J/ip decays with 3 billion ip(2S) sample, but the 
statistics will be substantially reduced. 

5 Charged-lepton flavor violating pro- 
cesses in decays of J/ip 

In the framework of the SM, lepton flavor is an 
accidentally conserved quantum number due to the 
massless neutrinos. The new physics beyond the SM 
introduces oscillations between the neutrino flavors 
by including the neutrino masses, which violate lep- 
ton flavor. However, SM processes involving charged- 
lepton flavor violation (CLFV) with massive neutri- 
nos are too tiny to be observed because they are sup - 



pressed by the quantity (AM 2 /M^) 2 < 10" 48 [41 1. 
Here AM„ is the difference between the masses of 
neutrinos of different flavor and M w is the mass of 
the charged weak vector boson. Hence, CLFV offers 
an opportunity to probe signature of new physics 42- 
The CLFV charmonium decays J/ip — !• IV (I 
and V = t, /i,e, I ^ V) are predicted by various the- 
oretical models that allow tree- level FCNC, includ- 
e.g., unparticle theory [46|, |47], i?-parity violat- 



ing, 



ing and large tan/3 SUSY scenarios, leptoquarks, and 
other models inspired by the idea of grand unifica- 

30- 



tion 



By using unitarity considerations, limits on CLFV 
jjL and r branching fractions [331 ] have been used to 



tions 



48 



■ 5fJ . In Ref. [48j , unitarity relations between 
the CLFV J/ip decay and the CLFV /i and r decays 
have been exploited. From the existing experimental 
bounds on the latter process [j^ . the stringent indi- 
rect limits can be obtained to be 4xl0~ 13 and 6xlfV 9 



for J/ip — > /ie and r/i/re, respectively. In Ref. [501 ] . 
the experimental bounds on nuclear /i~ — e~ conver- 
sion are used to place indirect limits on the CLFV 
J/ip decays, which are at the same level as those in 
Ref. [3] . As suggested by Zhang in Ref. [i^ , search- 
ing for the CLFV J/ip decays with a huge sample at 
the BESIII remains a worthwhile experimental chal- 
lenge. With a 58 M J/ip event sample at BESII, the 
following upper limits have been established 51| : 



B7l(J/ip -)• t^) < 8.3 x 10~ 6 ; 
BTZ(J/ip rV^) < 2-0 x 10" 6 ; 
BK(J/ip -» ^e^) < 1-1 x Ur 6 . 



(11) 
(12) 

(13) 



The limits on the two-body lepton flavor violating de- 
cays of the J/ip could be reduced to the 10 -8 ~ 10~ 9 
level at BESIII with one year full-luminosity run at 
the J/ip peak. This would be a significant improve- 
ment. 

We should note that the indirect limits on CLFV 
J/ip decays are significantly more stringent than the 
expected bounds from the BESIII experiment. How- 
ever, searches for the CLFV decays of vector mesons 
remain an important experimental effort since their 
observation at the rates above the indirect limits 
would be a manifestation of new CLFV physics, which 
does not fit into these theoretical analyses. In par- 
ticular, it may imply a nontrivial mechanism of self 
cancellation as discussed in Refs. 48 -5(3], which are 
considered as unnatural. 

6 Summary 

A mini-review of J/ip and ip(2S) rare decays in the 
BES-III experiment has been done. With one year's 
luminosity data, about 10 billion J/ip and 3 billion 
ip(2S) decays can be collected at BESIII experiment. 
The searches of weak decay, LFV decay, C or P vio- 
lated decay and invisible decays of charmonium states 
will be feasible. These rare decay processes will pro- 
vide opportunity for new physics space, especially, in 
charm sector. It should be noted that the recent mea- 
surement of 3.5cr evidence for direct CP violation in 
D° decay at LHCb may indicate possible contribu- 
tions to charm decays from physics beyond SM. We 
expect to see surprises from BESIII experiment in 
rare charmonium decays. 
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